Entanglement properties and momentum distributions of hard-core anyons on a ring. 
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We study the one-particle von Neumann entropy of a system of N hard-core anyons on a ring. 
The entropy is found to have a clear dependence on the anyonic parameter which characterizes the 
generalized fractional statistics described by the anyons. This confirms the entanglement is a valu- 
able quantity to investigate topological properties of quantum states. We derive the generalization 
to anyonic statistics of the Lenard formula for the one-particle density matrix of N hard-core bosons 
in the large N limit and extend our results by a numerical analysis of the entanglement entropy, 
providing additional insight into the problem under consideration. 



In recent years an intense research activity has been de- 
voted to the study of entanglement in many-body states. 
Initially, this effort has been mostly motivated by the 
fact that quantum correlated many-body states, which 
appear in various solid-state models, can be valuable re- 
sources for information processing and quantum compu- 
tation 0, Q ■ The theory of entanglement is now attract- 
ing even more attention because of its fundamental im- 
plication for the development of new efficient numerical 
methods for quantum systems [E 0, HI and for the char- 
acterization of quantum critical phases [E H 0] ■ 

Generally speaking, entanglement measures nonlocal 
properties of composite quantum systems and it can pro- 
vide additional information to that obtained by investi- 
gating local observables or traditional correlation func- 
tions. In this respect entanglement might be a sensi- 
tive probe into the topological properties of quantum 
states. A particularly significant quantity is the entan- 
glement entropy Sa, which is defined in a bipartite sys- 
tem A U B and quantified as the von Neumann entropy 
Sa = — TrpA In pa associated to the reduced density ma- 
trix pa of a subsystem A. In two-dimensional systems 
a firm connection between topological order and entan- 
glement entropy has been established in [E E3], where 
the entanglement entropy was defined by spatial parti- 
tioning. Recent studies on Laughlin states jll|, E2] have 
considered the entanglement entropy associated to par- 
ticle partitioning 11, EH- Also in this case, the entan- 
glement entropy turns out to reveal important aspects of 
the topological order in Laughlin States. 

The two-dimensional case is of particular interest due 
to the existence of models whose elementary excita- 
tions exhibit generalized fractional statistics. Anyons, 
the particles obeying such statistics, play a fundamental 
role in the description of the fractional quantum Hall 
effect [3]. Although this concept is essentially two- 
dimensional, anyons can also occur in one-dimensional 
(ID) systems [II EE EE El EE EE H3] , where statistics 



and interactions are inextricable, leading to strong short- 
range correlations. The ID anyonic models have proven 
useful to study persistent charge and magnetic currents 
in ID mesoscopic rings (Tsj . This possibility and their 
own pure theoretical interest lead us to investigate the 
effects of the anyonic statistics on the entanglement en- 
tropy in the present Letter. A discussion about quantum 
statistics and entanglement in a two-fermions system was 
introduced in 2l| and extended to the case of two bosons 
in (2^ . A mechanism of spin-space entanglement trans- 
fer based on the indistinguishability of two particle was 
proposed in [23| and shown to depend on the statistics 
(either fermionic or bosonic) of the particles involved. A 
comparison between the entanglement of two bosons and 
two fermions confined in a ID harmonic trap has been 
presented recently in [2^ |. 

In this Letter, we consider a system of N hard-core 
anyons on a ring which is the direct anyonic generaliza- 
tion of the Tonks-Girardeau gas. It offers a convenient 
framework to study topological effects: the many body 
ground state is known and its behaviour under the ex- 
change of two particles interpolates between bosons and 
fermions. We carry out an analytical and numerical anal- 
ysis of the dependence of the one-particle von Neumann 
entropy on the statistical parameter which determines 
the symmetry of the many body state. We derive the 
large N asymptotic expression of the anyonic one-particle 
density matrix. This asymptotic for m g eneralizes the 

25, 



and provides 



one obtained for hard-core bosons 
a one-parameter family of zero temperature momentum 
distributions interpolating between hard-core boson and 
free fermion distributions. Our results show that parti- 
cle entanglement depends in a non-trivial manner on the 
statistics and as such, may prove to be relevant to the 
study of topological properties of many-body quantum 
states. 

Let us consider a ID system of anyons confined on 
a ring of length L interacting each other via a repul- 
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sive 6— function potential. The model is denned by the 
Hamiltonian: 



H 



N 



dx 



2+7 E 

l<i<j<N 



(i) 



The TV— anyons wave function ty e (xi, x%, ..,xn) exhibits 
a generalized symmetry under the exchange of particles: 
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where e(x) = (or 1) if x > (x < 0) and is the 
anyonic parameter, defined as in Ref. [15j]. For = 
this model describes free fermions while, for 6 = n, it 
reduces to the Lieb-Liniger Bose gas. 

As first discussed in |l6[, the problem of ID anyons 
with contact interactions allows for an exact Bethe ansatz 
solution which shows that the Hamiltonian JT]) has the 
same energy spectrum than a ID interacting Bose gas 
with anyonic statistics-dependent effective coupling in 
the moving frame. Very recently, a detailed analysis of 
the low-e nergy p roperties of this model has been carried 
out in fl7L 1 1 8L ligj l . It was shown that the low-temperature 
thermodynamics of ID anyons with a <5-function potential 
coincides with the one of a gas of ideal particles obeying 
Haldane statistics: the interplay between the anyonic pa- 
rameter 9 and the coupling constant 7 determines a con- 
tinuous range of these generalized statistics. These stud- 
ies have shown that, for strong coupling, the dispersion 
relations of the anyon gas remain linear in the thermody- 
namic limit and the finite size corrections of the ground 
state energy found a central charge c = 1. 

In the case of spatial partitioning, the subsystem A be- 
ing a block of size I, conformal field theory results [27l.[28j 
predict the entanglement entropy (block entropy) Sa(1) 
to scale as Sa (0 ~ |lnL The dependence on the cou- 
pling constant and on the anyonic parameter is expected 
to show up in the sub-leading terms which, to the best of 
our knowledge, remain unknown. Below we will demon- 
strate that the one-particle entanglement entropy Sf (A) 
of A— anyons depends clearly on the anyonic parameter. 
Furthermore, we will show that, for a great number N of 
particles, the dependence of the entanglement entropy on 
the anyonic statistics appears in the sub-leading term. 

Let us consider the limit of hard-core anyons, i.e. 7 — > 
00. As recently shown in [2(|, the Fermi-Bose mapping 
method for one dimensional hard-core bosons [12] can be 
generalized to an anyon- fermion mapping (AF). Imposing 
the exclusion principle, i.e. the vanishing of the many- 
body wave function when two particles occupy the same 
position, the AF mapping reads [201 ] : 



* (xi, • • • ,xat) = 
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where ^q(xi, • • • , xn) is the N free fermion ground state 
function and A(xi,Xj) = e lSe ( Xi,x j\ In the following, we 
restrict ourselves to the case where N is odd, which cor- 
responds to a non-degenerate ground state. The topo- 
logical properties of the A-anyon wave function are en- 
coded in the factor Y\i <i< j <N A{xi,Xj) which gives the 
statistical phase e l9P resulting from the P exchanges 
needed for the particle positions to be brought to the 
ordering < Xi < X2 < • • • < xn < L. The peri- 
odic boundary conditions \& (xi, • • • , Xj + L, ■ ■ ■ , xn) = 

(xi, • • • , Xi, ■ ■ ■ , xn) impose the anyonic parameter to 
be an integer multiple of 2tt/ (N — 1), 6 = -^zjn. 

We are interested in computing the entropy Sf(N) = 
— Tr(/5jy lnpjy), where p%(x — x') is the one-particle re- 
duced density matrix: 



p e N (x - x') = 



l n 

Y\_dxi [* e (x,x 2 ,- • • ,xjv) 
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x* e (x',x 2 ,--- ,x N )\ , (4) 



normalized such that Pn{®) — !• I n momentum space, 
the entanglement entropy simply reads: 



S((N) 
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where c e N (n) = 1/L L p e N (x) cos(2n/Lnx) is the mo- 
mentum occupation in the ground state. For 8 = 
(free fermions), one has c%(n) = 1/N for (N — l)/2 < 
n < (N + l)/2 leading to the well known-result for free 
fermions = In A. The von Neumann entropy actually 
measures the uncertainty about the quantum state to at- 
tribute a state to the subsystem in consideration, which 
in this case comes exclusively from the fact the fermions 
are indistinguishable. 

The mapping §3§ between Fermi and anyon eigenfunc- 
tions preserves all scalar products and thus the energy 
spectrum and all the probability distributions involving 
the norm of the wave function are identical. Nevertheless, 
c e N (n) strongly depends on 9 as can be guessed from the 
drastic difference between the momentum distributions 
of free fermions and hard-core bosons. 

Analytical expressions of c 7r (n) for small (|n| -C 1) and 
large (|n| 3> 1) momenta have been found in the ther- 
modynamic limit [H, Hi], [35|, [3^, [I?! • These results show 
a In^ 1 / 2 singularity at n = 0, reflecting the tendency 
towards Bose-Einstein condensation. The corresponding 
result for finite N is more cumbersome. Using the N 3> 1 
asymptotic result for p 7T {x) [25], 



p 7T (x)^ Poo N 2|sin7rx/L|" 
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with poo = G(3/2) 4 /v / 2 and G(z) the Barnes G-function 
[40l | . c n {N) for N 3> n was shown to behave like [2^ |: 
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where T(z) is the standard Gamma function. We general- 
ized the above results to anyonic statistics. Representing 
the density matrix (0| in terms of a Toeplitz N— 1 x N— 1 
determinant: 



n Pn( x ) = detAr_i [<j) k j] (x) 



(8) 



where 
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we were able to compute the asymptotic form of p N using 
the Fisher-Hartwig conjecture 3a, 3t|. For N 1 the 
one-particle anyon density matrix reads: 
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and a{9) — ^(1 — We see that Eq. ([6]) is recovered 
for 9 = it. Details of the derivation and a more com- 
plete discussion of this result will be presented elsewhere 
4ll I42] ] . In the case of the generating function ((9j) the 
analysis of the behaviour of the Toeplitz determinant is 
subtle and the Fisher-Hartwig formula remains a conjec- 
ture. The validity of Eq. (fT0|) has thus be compared to 
the numerical evaluation of the determinant ([8} for finite 
N. For small 9, the convergence to the asymptotic result 
is quite slow and the formula provides a rough estimate 
for finite N, as can be seen for N=121 and 9 = ir/QO 
(Fig. [1] (a)). The likeness increases greatly with 9 and 
the agreement is already perfect with 9 — ir/2 for N = 61 
(Fig. [l](b)) and for 9 — 97r/10, close to hard-core bosons, 
for N=21 (Fig. ffl(c)). 

The Fourier coefficients of (TTU1) can be computed an- 
alytically [44|. The asymptotic behaviour of c n (N) for 
N n reads: 



e n (N) 
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where ri = n - [N(6/(2tt) - 1/2) + 0/(4tt) + 3/4J with 
\_x\ being the integer part of x. A comparison between 
the above formula and the exact results is shown in Fig. 
[21 Note that, from Eq. (fT0|) . the number of particles oc- 
cupying the low-energy state (n' = 0) scales as N"™' 1 , 
thus ruling out the possibility of anyon-condcnsation pre- 
dicted in the case of free anyons (4^|. For n N, 
the bosonic momentum distribution decays like n 
MS, S3- Based on exact results for N = 5,7, 9 and 



on numerical computations (see Fig. [2]), we expect this to 
be true for anyonic statistics as well. Under this assump- 
tion, the terms c 6 n (N) for n 3> N will not contribute 




FIG. 1: (Color online) Comparison between the numerical 
computation (red/dark gray curve) and asymptotic equation 
(black curves) for the one-particle anyon density for: (a): N = 
121,0 = 7r/60, (b): N = 61,6* = tt/2 and (c): N = 21, 9 = 
9tt/10. 



significantly to the entanglement entropy. While the 
main contribution to Sf(N) can be extracted from pip , 
the numerical simulations show that, for finite N, the 
crossover region between the asymptotic ([lip and power- 
law can not be neglected, especially near the fermionic 
point 9 = 0. 

We have determined numerically the value of the one- 
particle von Neumann entropy Sf(N) for different values 
of N (Fig. [3] (a)-(b)) which can be shown to scale with 
N in the following way: 



SZ(N)n\nN + f(0) + 



k(9) 



N 



(12) 



The anyonic-parameter-dependent function f(9) was de- 
termined numerically for system sizes for which the last 
term of this expansion can not be neglected. However, in 
the thermodynamic limit, only f(9) will remain relevant 
in our discussion. The result of our numerical analysis 
are displayed in Fig. [3[ where panel (a) is the data as 
fitted by Eq. ([T2")) for three values of the anyonic param 
eter and panel (b) is the resulting /(#). This function 
decreases monotonically from free fermions to hard-core 
bosons where it respectively takes the value /(0) = and 
s(tt) w —0.3, and from Fig. [3]-(b) it can be seen to be 
well fitted by a sinus function. 

To conclude, we have investigated analytically and nu- 
merically the one-particle von Neumann entropy and the 
momentum distributions of N hard-core anyons on a ring. 
We have determined the asymptotic expressions for the 
one-particle density matrix and for the momentum dis- 
tributions. Numerical results show the entanglement en- 
tropy exhibits a simple and non-trivial dependence on the 
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anyonic parameter even for TV — + oo, making it a suitable 
tool to study the topological properties of many-body 
quantum states. 
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FIG. 2: c n (N) obtained from numerical (Color online) and 
Eq.CIT} (Black line) versus n for = tv/2 and N = 61. The 
inset is the same plot for n > in a log-log scale (red/dark 
gray curve), while the dashed curve is a visual guide propor- 
tional to n -4 . 
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FIG. 3: (Color online) (a): Entanglement entropy as a func- 
tion of N for 6 = 7r (pluses), 9 = n/2 (crosses), 9 — n/10 
(dots) in log- linear scale, fitted according to Eq. (|12[1 . (b): 
f(9) obtained by numerical integration (pluses) and the cor- 
responding sinus fit (plain line). 
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